Universality classes of the Kardar-Parisi-Zhang equation 
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We re-examine mode-coupling theory for the Kardar-Parisi-Zhang (KPZ) equation in the strong 
coupling limit and show that there exists two branches of solutions. One branch (or universality 
class) only exists for dimensionalities d < dc — 2 and is similar to that found by a variety of analytic 
approaches, including replica symmetry breaking and Flory-Imry-Ma arguments. The second branch 
exists up to dc = 4 and gives values for the dynamical exponent z similar to those of numerical studies 
for d > 2. 
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The celebrated Kardar-Parisi-Zhang (KPZ) equation 
m was initially derived as a model to describe the ki- 
netic roughening of a growing interface and has been the 
subject of a great number of theoretical studies [3|- This 
is because the original growth problem has turned out 
to be equivalent to other important physical phenomena, 
such as the randomly stirred fluid (Burgers equation) 3| , 
directed polymers in random media (DPRM) |4|, dissi- 
pative transport p, |6[ or magnetic flux lines in supercon- 
ductors Q- The KPZ equation has thus emerged as one 
of the fundamental theoretical models for the study of 
universality classes in non-equilibrium scaling phenom- 
ena and phase transitions [^. 

It is a non-linear Langevin equation which describes 
the large-distance, long-time dynamics of the growth pro- 
cess specified by a single- valued height function /i(x, t) on 
a d-dimensional substrate x g IR : 

dth{x,t) ^iyV^h{x,t) + A/2 (V/i(x, t))^ + r]{x,t), (1) 

where r7(x, t) is a zero mean uncorrelated noise with vari- 
ance (?7(x,t)r/(x',i')) = 2DS'^{x-yi')Sit-t'). This equa- 
tion reflects the competition between the surface tension 
smoothing force vV^h, the preferential growth along the 
local normal to the surface represented by the non-linear 
term and the Langevin noise rj which tends to roughen 
the interface and mimics the stochastic nature of the 
growth. 

The stationary interface is characterized by the two- 
point correlation function C(|x — x'|, i — i') = ([/i(x, t) ~ 
/i(x',t')]^) and, in particular, its large-scale properties 
where C is expected to assume the scaling form C{d, r) — 
d^^ /(i?); where d = |x — x'| , r = \t — t'\ and x a-nd 
z are the roughness and dynamical exponents respec- 
tively. These two exponents are not independent since 
the Galilean symmetry [3| — the invariance of Eq. ([1]) 
under an infinitesimal tilting of the interface — enforces 
the scaling relation z + x = 2 for solutions associated 
with any fixed point at which A is non-zero. 

While some exact results are available in d = 1, yield- 
ing X = 1/2 and z = 3/2 0, 0, S], the complete theo- 
retical understanding of the KPZ equation in higher di- 
mensions is still lacking. For d > 2, there exists a phase 



transition between two different regimes, separated by a 
critical value Ac of the non- linear coefficient [l|, Isj. In 
the weak-coupling regime (A < Ac), the behavior is gov- 
erned by the A = fixed point — corresponding to the 
linear Edwards- Wilkinson equation [2] — with exponents 
X = {2 — d)/d and z = 2. In the strong-coupfing (rough) 
regime (A > Ac), the non-linearity becomes relevant and 
despite considerable efforts, the statistical properties of 
the strong-coupling (rough) regime for d > 2 remain con- 
troversial. 

The existence of a finite upper critical dimension dc at 
which the dynamical exponent z of the strong-coupling 
phase becomes 2 is also much debated. Most of the an- 
alytical approaches support a finite dc, but their pre- 
dictions are varied: one of the solutions of the mode- 
coupling (MC) equations and other arguments indicate 
dc — 4 [9|, [lOJ, [ill, ll3|, whereas functional renormaliza- 
tion group to two-loop order suggests dc — 2.5 [l3|. A 
set of related theories such as replica symmetry break- 
ing [ij], variational studies 15 1 and Flory-Imry-Ma ar- 
guments [l6| give dc = 2. On the other hand, numerical 
simulations and real space calculations find no evidence 
at all for a finite 4 [llHIisl! 



We re-examine in this paper MC theory. It is a self- 
consistent approximation [5|, |6|, where in the diagram- 
matic expansion for the correlation and response func- 
tions only diagrams which do not rcnormalize the three- 
point vertex A are retained. The MC approximation has 
been widely studied 0, B [ul l20i |2l| . Furthermore the 
MC equations are exact for the large A^-limit of a gener- 
alised iV-component KPZ [20|. This allows in principle 
for a systematic expansion in 1/A^. We shall show that 
within the MC approximation there exists two solutions, 
one of which had been previously overlooked. The new 
solution (universality class) only exists when d < 2 and 
seems similar to the solution found in the analytical stud- 
ies which give dc = 2 [3, [H, |ll| . The other MC solution 
exists on the whole range < d < 4 and has dc = 4. 
This previously known solution gives in d = 2 a value for 
z close to that from numerical work, e.g. [18,] . 

The correlation and response functions are defined 
in Fourier space by C{k,uj) ~ {h{k,uj)h*{k,uj)) and 



G{k,oj)5''{k~k')5{u^^uj') = (lll^l^), where (•) in- 

dicates an average over 77. In the MC approximation, the 
correlation and response functions are the solutions of 
two coupled equations, 



G-'(.,.)^G.-.(.,.) + A^/f/(0, 



q- (k-q) 



C{k,uj) 



X {q-k)G{k~q,UJ~-n)C{q,n) (2) 

= C,{k,.) + -\G{Ku)\ J ^J j^ 
X {q-{k~q)fC{k-q,oj-n)C{q,n) (3) 



where Gf){k,uj) = (vk^ — iu) ^ is the bare response 
function, and Co(fc,w) = 2D|G(fc,a;)p [1]. In the scal- 
ing hmit, G{k,uj) and C{k,uj) take the forms G{k,uj) — 
k-'^g (uj/k'') and C(fc, w) = fc-(2x-H<i-H^)ji (u;/fc2) and Eqs. 
([2]) and ([3]) translate into coupled equations for the scal- 
ing functions n(x) and g{x): 



n{x) 



—ix + Ii{x) , 
\9i^)\'l2{x), 



(4) 
(5) 



where x = co/k^ and /i(3:) and l2{x) are given by [9| 

/•TT /'OO 

/i(a;) = ^/ dOsin'^-^e dq cos 6 (cos 9 - q) 
Jo Jo 

y^q'^-'r-^ p dyg(^—pL^n{y), 

-^2(2;) = -/ dOsin'^-^e I dq{cos0^qf 
2 Jo Jo 

with P = AV(2'^r(^)7r(''+3)/2)^ ^2 ^ l + q2_2qcos6'. 
Notice that the "bare term" Co{k,uj) has been dropped 
from these equations as it is negligible in the scaling limit 
provided 2x + d + z > 2z, i.e. 4 + d > 3z. 

All the (necessarily approximate) solutions of the MC 
equations involve an ansatz on the form of the scaling 
functions n{x) and g{x). The relation z = z{d) is then 
obtained requiring consistency of Eqs. (|3]) and © on 
matching both sides at an arbitrarily chosen value of x. 
Due to the non- locality of Eqs. @ and (O, the match- 
ing condition depends on the form of the functions n and 
g for all a;, so the ansatz needs to be reliable for all x, 
and in particular satisfy the large x asymptotic forms 
n{x) ~ x"^^''/^, gB.{x) ^ x"^^"^!^ and gi{x) —> x^^ 
where (3 — d + 4 — 2z and g{x) — gnix) + igi{x). Co- 



laiori and Moore [ll|, II4I proposed an ansatz which sat- 
isfies these large x constraints. It enabled the authors 
to provide numerical estimates of z when d = 2, 3 in rea- 
sonable a gre ement with exponents obtained from simula- 
tions 171, [l8| . It also yielded an integer finite upper crit- 
ical dimension dc = 4 lll|. Moreover, it led in d = 1 to 



the discovery of a stretched exponential behavior for the 
two-point correlation function 1^ similar to that found 
in the exact solution |22l |. 

Using a similar ansatz, we next show that Eqs. @ 
and (O have an additional solution with dc = 2, before 
discussing its significance. The advocated ansatz is most 
conveniently expressed in Fourier space, 

e{p)eM-\p?''-), 



dip) 

n{p) 



= Aexp{~\Bpf^") 



(6) 

(7) 



where n and g are the Fourier transform of n{x) and 5(2;) 
respectively, A is a (nonzero) parameter and B = B(d^ z) 



an arbitrary function 11|. Requiring the matching of 



Eqs. (dl and ([5]) at large x then yields: 

GdAB-^ = d{2-z)/{z^I{B,d,z)) 
GdAB-^ = {l3lzf{p-z)2'^/TC'^ 



P 



(8) 
(9) 



where Gd = \'^T{{d - l)/2)/(4 7r'^ r(d - 1)) and 

/•OC 

I{B, d,z)^ dsB{l~ 2s^)s^'-^ exp ( - B^'/'sf' - s^) . 

_ , Jo 
Combining ([8]) and ([9]) then gives an implicit equation 

for the curve z(d) depending on B only: 



/?' (/? - z) 



2(2z-0)/P 



HB,d,z) 



d{2- z) 
B 



(10) 



Eq. ((T0|) can be solved numerically once B{d, z) has been 
fixed. Prior to discussing the choice of B, we stress the 
(_B-independent) intrinsic properties of the solutions Sd.z- 
First, as z -^ 2, I{B, d, 2) is finite for any _B > 0. One 
can then infer from Eq. pop that 82,2 and ^4.2 are al- 
ways solutions independent of i? (provided B~^ is finite). 
This singles out the two dimensions dc = 2 and dc = 4 
as critical ones since z — > 2 for d -^ d^. Second, as 
d -^ 0, I{B, 0, z) is non-zero for any B > 0, and Eq. 
pop yields the two solutions S'0,4/3 and S'0,1 (provided 
again B~^ is finite). Lastly, as the MC equations satisfy 
the fluctuation-dissipation theorem (FDT) that exists in 
d — I [2l| yielding the exact result z = 3/2, we expect 
the ansatz to preserve this property, which is fulfilled 
for 5(1,3/2) = 1 as can be checked on Eqs. ^ and 
([7]). It follows that S*! 3/2 is solution of pop . To summa- 
rize, Eq. pop has five intrinsic solutions Sd,z which do 
not depend on the choice of B (provided B~^ < 00 and 
5(1,3/2) = 1). 

Furthermore, Eq. ([5]) provides additional constraints 
on B{d, z) in the two limits z ^ 2 and d ^ 0. As 
I{B, d, 2) is finite as z — > 2, the right-hand side of ([8]) 
vanishes identically, which forces B~^ ~ (2 — z) on the 
left-hand side. One can then derive the expansions of two 
solutions, labelled zs and z-p, about their critical dimen- 
sions: 

ZF = 2-(4-d)/4 + 0((4-d)2) (11) 

zs = 2-(2-d)/3 + 0((2-d)2). (12) 



The arguments of Ref. llj can be used to show that 
these expansions are very robust and more general than 
the particular ansatz of Eqs. ^ and ([7]). Similarly, as 
I{B,0,z) is finite for d — ^ and z > 1, Eq. ([5]) imposes 
B~^ ~ d. The expansions of the two solutions about 
d — can also be worked out: zp — 4/3 + d/3 + 0{d^) 
andzs = 1 + ^/4 + 0(^2). 

We emphasize that all these features do not depend 
on the choice of B. Any ansatz would lead to the same 
qualitative picture as Fig. [T]provided B ^ {d{2 — z))~^ . 
The specific one used here further exploits the uniqueness 
of the solution in d = 1 which imposes additional con- 
straints on B [26|. With this choice, we have solved Eq. 
(fTU)l numerically in < d < 4 and found that there exists 
two branches of solutions za,[d) and z-p{d) lying in the in- 
terval d e (0, 2) but only one solution, the continuation 
of zp, in the interval d € (2,4). 




FIG. 1: The two solutions zs(d) and 2p(d) of the MC equa- 
tions found in this work, compared with numerical 18] results 
for d = 2, 3, 4 and the exact results in d = 0, 1, all represented 
by squares. 

The branch zp, referred to as the F solution as it exists 
over the full dimensionality range d £ (0, 4), is similar to 
the solution analyzed in I 111 27 1. It lies fairly close to 
typical numerical results [l8| with z — 1.632 for d = 2 
and z = 1.787 for d = 3, reproduces the exact results in 
d = and d = 1 (see Fig. [T]) and reaches z = 2 at dc = 4, 
which defines the upper critical dimension of this solu- 
tion. The other branch, denoted S as it is associated with 
the shorter dimensionality interval d £ (0,2), embodies 
a new universality class with zs ^ zp for d ^ 1. For this 
solution, z = 2 at dc = 2 and it does not exist for d > 2, 
so that there remains a single strong-coupling solution in 
and above two dimensions - the F solution. We enrpha- 
sise that both the F and S solutions are non-perturbative 
since no (stable) perturbative fixed point can be found 
for d < 2 within the standard perturbative expansion in 
the coupHng constant g = 2}?D/u^, even to all orders 
[23'|. MC theory involves resumming an infinite set of 
diagrams and so is capable of going beyond standard RG 



techniques. 

We now comment on our findings. According to MC 
theory, for d < 2 there are two possible strong-coupling 
universality classes, F and S. It is clear that for a given set 
of initial conditions and bare correlators, the time evolu- 
tion would unambiguously steer the system towards one 
or other of the two solutions. Equivalently, in the map- 
ping to the DPRM, the solution with lowest free energy 
would be selected as the physical one. The apparent free- 
dom to choose a solution arises because we have consid- 
ered the scaling regime in which the bare correlators have 
been dropped as irrelevant, so we have lost track of their 
role. However, a clue as to the form of bare correlator 
needed for the S solution to emerge can be obtained by in- 
specting the zero-dimensional case. In d = 0, the DPRM 
analogue is a directed walk along a chain of length L on 
which the site energies are random. A dynamical expo- 
nent z = 1 would correspond to free energy fluctuations 
of order L^/z-i _ ^ which is larger than the \/Z (for 
z = 4/3) expected for independent random site energies, 
with z = 1 could only exist with long-range correlations 
between the site energies and is therefore not the natu- 
ral solution for the short-range problem in d = 0. We 
suspect that the zs branch in the whole range < d < 2 
would only be realized if there were some long-range cor- 
relations in the noise for the KPZ problem or between 
the site energies for the DPRM analogue. 

To investigate this possibility further, suppose instead 
of starting out with noise with short-range correlations, 
one introduces some long-ranged correlated noise, of the 
form {rj{x,t)rj{x\t')) = 2R{x — x')6{t — t') with some 
power law distribution R{x) oc x'^p~'^, or equivalently in 
Fourier space R{k) = D{1 + wk^'^P). The MC equa- 
tions would formally be the same with D replaced by 
the 'generalized' noise correlator Rik) [2l|. These long- 
range bare noise correlations will not affect the scal- 
ing behavior of the system as long as they decay suf- 
ficiently quickly, so that they become irrelevant after 
some transient regime. The condition for this to hap- 
pen can be worked out by substituting the scaling forms 
for G and C into Eq. ^. The initial correlations de- 
cay as Co(fc,w) = 2i?(fc)|G(fc,w)|2 - fc-2(2=+p) whereas 
the long-time correlation function is expected to behave 
as C(fc,w) ~ fc-(2x-(-'i+2=) = }^-(d+i-z) ygj^g ^jjg identity 

z + X = 2. Comparing these two expressions, one infers 
that the long-range correlations will destabilize a short- 
range solution with dynamical exponent z when 



P> 



3z 



(13) 



which is identical to the criterion obtained in 24] from 
RG arguments. The condition of Eq. (|13p determines 
the phase boundary p(d, z) between the short-range and 
long-range stability domains for each of the two MC so- 
lutions Zp and zg. In the long-range noise dominated 
region zlr = <>-+i-2p 24] xhe solutions F and S are sta- 



ble against noise of range p if their associated z satisfies 
z < 2lr- For all < d < 2, we believe that the S solution 
requires for its existence the addition of long-range noise 
of infinitesimal amplitude with pg = '^^'^"'^^'^ , as in the 
case d — (where z = 1 would correspond to long-range 
noise with p = 0.5^). For < d < 1, the S solution 
lies at the stability limit of the long-range solution and 
in the interval 1 < d < 2, it will be the solution to the 
KPZ problem in the presence of infinitesimal amounts of 
such correlated noise if (say) the associated free energy in 
the DPRM analogue is lower than that of the F solution 
(with the same correlated noise). This S solution could 
hence possibly be brought out by tuning the parameters. 
We believe that in the absence of some long-range corre- 
lation, the F solution will be the appropriate solution or 
universality class for < d < 4. 

Our S solution might appear to be of academic inter- 
est as it only exists when d < 2 and coincides with the 
F one in d = 1. However, it turns out to be strikingly 
similar to the solution found in several theoretical stud- 
ies. We focus on three of them: the replica symmetry 
breaking (RSB) calculation of Mezard and Parisi [14|, a 
Gaussian variational calculation ^5|] and Flory-Imry-Ma 
(FIM) arguments [16j. All three approaches result in a 
strong-coupling solution with dc = 2. The RSB treat- 
ment yields a dynamical exponent equivalent to the FIM 
result [il: zfim = (4-f-d)/3 ~ 2-(2-d)/3 + C'((2-d)2), 
which is identical to (|12p and suggests that the MC so- 
lution S and the RSB - FIM solution are indeed differ- 
ent approximations to the same underlying universality 
class. It is not apparent how the long-range correlations 
needed for the S solution to exist enter these analytical 
approaches. However, it has been observed long ago that 
by the addition of a type of long-range force, the Flory 
formula in the conventional self-avoiding polymer prob- 
lem becomes exact ^2S] . 

To summarize, we have found a new MC solution, S, 
which has an upper critical dimension dc = 2, in addition 
to the usual one, F, which has dc ~ 4. This new solution 
seems similar to that found by analytical treatments such 
as RSB and FIM arguments, but we believe will only 
exist in the presence of appropriate long-range correlated 
noise. We acknowledge that the numerical value for z 
shown in Fig. [ijwhen d = 4 is not apparently consistent 
with the upper critical dimension dc — 4 of the F solution, 
but we attribute that to the inability of numerical studies 
in high dimensions to access the scaling regime. 
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